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We show that, for fermionic atoms in a one-dimensional optical lattice, the fraction of atoms in 
doubly occupied sites is a highly non-monotonic function of temperature. We demonstrate that 
this property persists even in the presence of realistic harmonic confinement, and that it leads to a 
suppression of entropy at intermediate temperatures that offers a route to adiabatic cooling. Our 
interpretation of the suppression is that such intermediate temperatures are simultaneously too 
high for quantum coherence and too low for significant thermal excitation of double occupancy thus 
offering a clear indicator of the onset of quantum fluctuations. 
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I. INTRODUCTION 



The international effort to emulate the behavior of cor- 
related electrons in solids using ultracold atomic systems 
^{5] is by now a familiar topic. Experimental progress 
has been rapid, and the Mott insulating state of fermionic 
atoms in a 3D cubic optical lattice has been observed 
[HIS]. However, to explore the key open questions about 
the Hubbard model — for example, how similar its phase 
diagram really is to that of the high-temperature super- 
conductors [5] — it is necessary to achieve a significantly 
lower entropy per particle than at present. To be quanti- 
tative, the entropy per particle must drop below fcs ln2; 
this is the scale at which antiferromagnetic correlations 
begin to set in. 

The role of theory in recent work on this problem has 
been vital. On the one hand, experimental work is being 
conducted to help elucidate the properties of the Hub- 
bard model at low temperatures. On the other, the inter- 
pretation of the data taken in experiments often depends 
heavily on the theoretical understanding of the Hubbard 
model itself. For example, one experiment [5] measured 
the density profile of the cloud under varying amounts of 
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harmonic confinement; the Mott insulator was identified 
by close comparison with theoretical predictions of the 
same. In more recent experiments [7], measurements of 
the fraction of doubly occupied sites in the optical lat- 
tice, d, were compared with theory to show that these 
experiments are still dominated by thermal (rather than 
quantum) fluctuations, a topic we take up below. 

As well as being relevant for interpretation of mea- 
surements on ultracold atomic systems, good theoretical 
understanding is crucial for the development of experi- 
mental cooling protocols. Here also the behavior of the 
double-occupancy fraction as a function of temperature is 
important, as it has a direct relation to the entropic prop- 
erties of the system. For example, the Pomeranchuk-like 
cooling method outlined in |8j relies on an entropy en- 
hancement, observable as a surprising suppression of dou- 
ble occupancy at intermediate temperatures. Whether 
this effect is strong enough to be useful is thus a key 
question; recent numerical calculations tackling this is- 
sue in 3D lattices [S] are suggestive but results in 2D 
lattices [101 E] are mixed. Highly accurate theoretical 
methods will be needed at low temperatures. 

Motivated by these requirements — interpretation of 
data and design of cooling protocols — we address in 
this work the following two questions. First, how good 
is the double occupancy as a probe of the temperature? 



And second, under what circumstances do we see the 
suppression of double occupancy at intermediate tem- 
peratures required for the adiabatic cooling protocol of 
Ref. [5]? We construct an accurate theoretical framework 
for treating low temperature fermions in one dimension. 
We shall find that, in a ID optical lattice, d{T) is a non- 
monotonic function of temperature. From the point of 
view of temperature measurement, this is a disadvan- 
tage; on the other hand, the resulting double-occupancy 
suppression indicates a large entropy enhancement, of- 
fering a clear route to adiabatic cooling. We shall argue 
that the enhancement of double occupancy occurs for 
different reasons in different temperature regimes: it is 
driven by thermal fluctuations at high temperatures, but 
by quantum fluctuations at low temperatures. We thus 
interpret the suppression in double occupancy at inter- 
mediate temperatures as occurring because the system is 
— loosely speaking — too hot for quantum coherence, 
but too cold for significant thermal excitation of double 
occupancy. Finally, we shall show that these local effects 
can be measured quantitatively with a bulk observable 
even in the presence of realistic harmonic confinement. 

Our method combines zero-temperature Bethe Ansatz 
(BA) studies with finite-temperature series expansions. 
Recent implementations of BA [15] compute observables 
such as the core compressibility i^6|. Here, however we 
show that double occupancy alone indicates a large en- 
tropy enhancement, even though the ID system is not 
a Fermi liquid. By focusing on ID Hubbard physics we 
demonstrate that exact calculations over the entire tem- 
perature range can be used as a platform to guide exper- 
iments in the construction of low temperature quantum 
states in optical lattice emulators. 
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Figure 1. Double-occupancy fraction versus temperature for 
the uniform ID Hubbard model at half filling for several U/t. 
The dashed line indicates the low-T expansion Kt; the solid 
line is the fourth-order high-T expansion; the dot-dashed line 
is the exact result for a two-site Hubbard model in the grand 
canonical ensemble; and the dotted line is the same but with 
the number of particles set to exactly two. At low T, non-zero 
d arises from coherent quantum fluctuations, but at interme- 
diate temperatures this coherence is destroyed, leading to a 
suppression of d{T) shown by the double-headed arrow. 



II. MODEL 

We consider the Hubbard model of cold fermionic 
atoms in a one-dimensional optical lattice: 



trapping, it would be given by 
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The indices i and j are integers labeling the sites of the 
lattice; the index a ='f,i labels two different hyperfine 
states of the A'^ atoms in question, and we take both of 
these states to have equal total populations N^ = N^ = 
N/2. The operator Cj^ annihilates a fermion of species a 
on site j; the number operator rij^ is defined as Cj^Cja- 
The first term in (II]) represents the quantum tunnel- 
ing of atoms between neighboring sites of the lattice; 
the second the repulsion due to s-wave scattering when 
two atoms are present on the same site; and the third 
the single-atom site energy. This last term may be site- 
dependent; for example, if the system were confined to 
a 'box' of length 2L we would have e^ = for ja < L, 
e, = oo otherwise. In the more realistic case of harmonic 



ta'f/L', 



(2) 



where a is the lattice spacing and L is a length scale 
related to the curvature of the trap. The hopping pa- 
rameter t depends on the depth of the optical lattice; 
the on-site repulsion U on the s-wave scattering length, 
which is tunable via Feshbach- resonance methods; and 
the on-site energies Ej on the shape and amplitude of the 
trapping potential. 



III. RESULTS (HOMOGENEOUS CASE) 

We begin our analysis with an analytic study of the 
zero temperature double occupancy in the bulk limit (ob- 
tained by assuming a 'box' potential as defined above and 
taking the thermodynamic limit L, A — ^ oo, N/L = con- 
stant) at half filling (/x = U/2). The double occupancy 
fraction, defined as the fraction of atoms in doubly occu- 



pied sites, is given by: 



d = {2 IN) {dF/dU)^^^ = {2 IN) J2 (njnj) , (3) 

i 

where F = E — TS is the free energy for entropy S and 
N is the total number of particles. We compute the bulk 
double occupancy fraction at half filling from the exact 
solution of Lieb and Wu [13] : 
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dujJo{uj)Ji{uj)sech'^{ujU/4:). (4) 



Here, J„ is the n*'^-order Bessel function. Here and in the 
following we work in units of t unless otherwise noted. 
This expression shows that strong quantum fluctuations 
in the ground state induce finite double occupancy, even 
at (n) — 1 and large U. 

Weak thermal fluctuations compete with quantum 
fluctuations at low temperatures (T <C t). Using results 
from the quantum transfer matrix method |17| we include 
thermal fluctuations to find the exact expression for the 
low T double occupancy fraction: 




Figure 2. Schematic of trap filling in three different regimes 
with decreasing U from left to right. Left: Mott insulator 
with single occupancy, Center: Doubly occupied sites at the 
trap center flanked by edges with low filling, and Right: A 
band insulator dominated by doubly occupied sites. These 
sketches ignore both quantum and thermal fiuctuations. 



d{T) > dujviU) ~ C{2n/U)T^ + OiT^^), (5) 
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where the function: 
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is related to the unity central charge predicted by confor- 
mal field theory for the Heisenberg universality class [17] . 
/ is the n"^-order modified Bessel function. It is striking 
to note that C > for all U. Thus weak thermal fluctu- 
ations counterintuitively lower the double occupancy in 
the ID Fermi Hubbard model. In contrast, large ther- 
mal fluctuations dramatically increase double occupancy. 

We use a fourth order high temperature series expan- 
sion of the ID Hubbard model to include large thermal 
fluctuations. Expansion of the free energy in powers of 
t/T about the atomic limit allows calculation of high 
temperature observables that compare well with finite 
temperature BA calculations [IH]. We compute the dou- 
ble occupancy fraction to fourth order. For brevity, we 
present the equation only up to second order (for the full 
expression see the Appendix): 
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Figure 1 matches the low T expansion (Eq. [s]) with 
the fourth order high T expansion of d within their 
respective regimes of convergence. Pronounced dips 
form when quantum and thermal fluctuations compete 
to dramatically enhance the entropy (The entropy can 



be obtained from the Maxwell relation {dS/dU)j, ^ = 
-{N/2) {dd/dT)jjj^). The finite value of d at T = is 
due entirely to quantum fluctuations but the double ar- 
row indicates the effect of thermal fluctuations at T = 2t. 
The related Pomeranchuk-like effect in higher dimen- 
sional lattices has been discussed in terms of Fermi liquid 
properties [HI [19] but it is well known that ID Hubbard 
models exhibit non-Fermi liquid behavior. 

We can understand the competition between quantum 
and thermal fluctuations in terms of a two site Hubbard 
model, which does not rely on Fermi liquid effects. At 
half filling the two particle sector has the lowest energy 
at T = for C/ > 2t. The ground state uses quan- 
tum fluctuations to lower its energy by hybridizing basis 
states. The singly occupied basis states correspond to 
the singlet (| t4) - I iA))/V2 and the triplet {| t,t), 
(I t,l> + I l,t>)/V2, I ;,;)} states. The singlet wiU 
mix with the doubly occupied states {| tijO): |0,ti)} to 
lower its kinetic energy without paying a large penalty 
from the interaction e nergy. The ground-state energy 
becomes ^o = U/2 - ./iP/A + W « -At^ /U for t < U. 

At small but finite temperature, the three triplet 
states, all with energy Ei — 0, allow thermal fluctua- 
tions to increase the entropy considerably (and therefore 
lower F) as we increase T. However, these 3 excited 
states involve no double occupancy and therefore d de- 
creases. As T is increased further, however, even higher- 
energy states are populated. These other states, e.g., 
I tiiO) — |0,ti), involve double occupancy again, so the 
double occupancy goes up once more for large T . The 
non-monotonic behavior of double occupancy is closely- 
related to the separation of energy scales between the 
spin and charge degrees of freedom [^TJ [55] . 
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Figure 3. The circles indicate double occupancy fraction ver- 
sus particle number in a trapped system (L — 50) at T = 
for half filling at the trap center and U/t — 8 computed using 
BA-LDA. The solid line is the bulk value in the Mott insulat- 
ing ground state (Eq. |4|. The insets, from left to right, show 
the corresponding density as function of position in the trap 
for 100, 216, and 300 particles, respectively. 



Fig. [T] plots d for the two-site Hubbard model for 
comparison with the thermodynamic limit. The pro- 
nounced dip in d corresponds to a finite T entropy gain. 
The remarkably close agreement between the expansions 
and the exact two-site results suggests that spinous and 
holons could be thought of in a localized picture. The 
close agreement also suggests that entropy gain when 
increasing interaction at intermediate temperatures is a 
generic feature of few-particle bonding in non-frustrated 
Hubbard models. 



IV. RESULTS (HARMONICALLY TRAPPED 

CASE) 

We now address the competition between thermal and 
quantum fluctuations under realistic parabolic trapping. 
Let us first consider the situation in the absence of fluctu- 
ations, sketched in Fig. [5] The left and right schematics 
show two scenarios where double occupancy should show 
a clear indication of the underlying trapped state. For 
large U (left) the deep Mott insulator regime corresponds 
to strong suppression of all doubly occupied sites. The 
low U regime (right) shows the band insulator, identifi- 
able with a large number of doubly occupied sites. The 
central schematic shows a mixture of both singly occu- 
pied sites at the edges and doubly occupied sites near 
the trap center. When fluctuations are turned on, the 
situation will be complicated further by the appearance 
of metallic regions and finite double occupancy at sites 
where (n) < 2. 

At first it may appear that trapping will hamper ef- 
forts to distinguish between thermal and quantum fluc- 
tuations in observations of double occupancy. We will 
use the Bethe-Ansatz local-density approximation (BA- 
LDA) [33] to show that bulk values are indeed observable 
in strongly inhomogeneous and relatively small trapped 
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Figure 4. The circles and solid line show the same as Fig. [3] 
but for three values of U. The dashed (dotted) lines show re- 
sults computed from an exact second order high temperature 
series expansion a,t T — 2t {T — 4t). The plateau at finite 
d indicates strong quantum fluctuations. The double arrow 
in the bottom panel indicates the observable effects of weak 
thermal fluctuations in a trap and corresponds to the double 
Fig.0 
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systems. Interestingly, local measurements [20] are not 
required to recover the bulk physics: a global measure- 
ment carried out on the whole trapped system suffices. 

We consider finite values of L in Eq. (21) to study 
trapped systems with up to iV ^ 500 particles with BA- 
LDA. Fig. [3] plots the double occupancy fraction versus 
Na/2L, an effective filling in the trap in the thermody- 
namic limit L,N ^ cx), Na/2L = constant [121 [13]. For 
low particle numbers, Na/2L < 1, the trap is entirely 
compressible due to edge effects. Near Na/2L ^ 2 the 
central region forms a Mott insulator at the trap center. 
At T = the system is dominated by quantum fluctua- 
tions manifest in the finite double occupancy even in a 
regime with (n^) < 1 for all i. d in the trap converges to 
duw in the formation of a plateau. Observation of this 
plateau indicates a strong Mott insulator with quantum 
fluctuations, d excludes edge effects and therefore allows 
a measurement of the bulk value of d, du/Vy even in an 
in-homogenous mesoscopic system. Increasing N further 
turns on double occupancy at the trap center leading to 
a pronounced cusp near Na/2L ^2.5. 

We use a 2"^^ order high temperature series expansion 
to include finite temperatures in the trap. The series 
from Ref. [H] is adapted to ID. The series we use for the 
trapped system is exact up to 0{{t/T)^). Fig.Hcompares 
the second order series at finite T and the BA-LDA at 
T = 0. From the top panel we see that strong thermal 
fluctuations tend to increase d above the bulk value for 
the Mott plateau near Na/2L ~ 2.3. For larger U the 
Mott gap suppresses thermal fluctuations. In the bottom 
panels we see that thermal fluctuations tend to decrease 



d in the trap for T < 2t. The double arrow Hne in the 
bottom panel shows that the effects of weak thermal fluc- 
tuations for T = 2t in a uniform system (Fig. IT]) are also 
observable in a trap. We predict that even lower tem- 
peratures will tend to increase d to its bulk value, c^lw 
in the Mott regime. Thus an increase of d with decreas- 
ing T demonstrates an observable capable of pinpointing 
a regime with dominant quantum fluctuations and low 
entropy per particle. This regime can be used for adia- 
batic cooling and constructing higher dimensional optical 
lattices with low entropy. 



CONNECTION TO ADIABATIC COOLING 



Our results complement a protocol for adiabatic cool- 
ing [5]. A Maxwell relation can be used to show that 
the thermal suppression of double occupancy implies a 
suppression of temperature with U [8]. Adiabatic cool- 
ing proposes to use changes in a tunable optical lattice 
parameter, U , to cross phase boundaries at fixed entropy. 
The suppression of double occupancy that we find here 
implies that one dimensional systems offer a controlled 
platform for adiabatic cooling that can be used to sys- 
tematically prepare low temperature systems in higher 
dimensions. Lowering the optical lattice depth along one 
(or two) perpendicular directions in the lattice adds the 
term —t± ^^ ,^ ■■. cl^c,^ + H.c. to the Hamiltonian, in- 
creasing the dimensionality from ID to 2D or 3D. Thus 
by isentropically changing t±, experiments will be able 
to identify and prepare higher dimensional optical lat- 
tice emulators from benchmark ID configurations at en- 
tropies per particle where quantum fluctuations domi- 
nate, i.e., below ks^n^. Although this does not put the 
system in a regime of low entropy per particle a priori, 
adiabatic cooling may allow a 3D system to cross the 
transition line from a paramagnetic to an antiferromag- 
netic phase. 



VI. CONCLUSION 

We have studied the behavior of the fraction of atoms 
in doubly occupied sites, d, for the homogeneous and 
trapped Hubbard model. In the homogeneous system, 
d shows a dip as function of temperature which signals 
the onset of antiferromagnetic correlations. We also find 
that this important regime can be clearly identified even 
in harmonically trapped systems. The double occupancy 
fraction increases with the number of particles in the 
trap, but remains fiat when the trap has a large central 
plateau in its density profile. The temperature behavior 
of d in the harmonic trap with a large central plateau 
is similar to that of half-filled homogeneous systems and 
can be used to conclusively identify a regime with domi- 
nant quantum fluctuations and low entropy per particle. 
Such an identification will set the stage for the prepara- 
tion of states with antiferromagnetic correlations in op- 
tical lattices. 

In preparation of this manuscript we recently became 
aware of related work in several different lattice geome- 
tries, Ref. [24j, that also supports a similar entropy en- 
hancement scenario to that put forward in Ref. [IQia but 
beyond the DMFT approximation. 
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VIII. APPENDIX 

Here we present an expression for the double occu- 
pancy fraction at half filling in units of i = 1 up to fourth 
order in t/T based on the free energy of Ref. [18] : 
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